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Abstract. Chaotic tunneling in a driven double-well system is investigated in ab- 
O ' sence as well as in the presence of dissipation. As the constitutive mechanism of chaos- 

l) , assisted tunneling, we focus on the dynamics in the vicinity of three-level crossings in 

the quasienergy spectrum. They are formed when a tunnel doublet, located on a pair 
of symmetry-related tori in the classical phase space, approaches a chaotic singlet in 
energy. The coherent quantum dynamics near the crossing, in particular the enhanced 
tunneling that involves the chaotic singlet state as a "step stone", is described sat- 
►^ \ isfactorily by a three-state model. It fails, however, for the corresponding dissipative 

(— ^ ■ dynamics, because incoherent transitions due to the interaction with the environment 

(^ ' indirectly couple the three states in the crossing to the remaining quasienergy states. 

^N) , We model dissipation by coupling the double well, the driving included, to a heat bath. 

c7^ ' The time dependence of the central system, with a quasienergy spectrum containing 

exponentially small tunnel splittings, requires special considerations when applying 

O^ , the Born-Markov and rotating-wave approximations to derive a master equation for 

the density operator. We discuss the effect of decoherence on the now transient chaos- 
assisted tunneling: While decoherence is accelerated practically independent of temper- 
ature near the center of the crossing, it can be stabilzed with increasing temperature at 
a chaotic-state induced exact crossing of the ground-state quasienergies. Moreover the 
asymptotic amount of coherence left within the vicinity of the crossing is enhanced if 
the temperature is below the splitting of the avoided crossing; but becomes diminished 
Q ' when temperature raises above the splitting (chaos-induced coherence or incoherence, 

^ , respectively). The asymptotic state of the driven dissipative quantum dynamics par- 

tially resembles the, possibly strange, attractor of the corresponding damped driven 
classical dynamics, but also exhibits characteristic quantum effects. 









X 



C^ ' 1 Introduction 



The interplay of classical chaos and dissipation in a quantum system bears inter- 
esting effects at the border between classical and quantum mechanics like, e.g., 
the suppression of classical chaos by quantum interference [1] or its restauration 
by dissipation [2] . While the mutual influence of quantum coherence and classical 
chaos is under investigation since many years, the additional effects caused by 
coupling the chaotic system to an environment, namely dissipation and decoher- 
ence, have been studied only rarely. One reason is that by including dissipation, 
the computational effort grows drastically, since one has to deal with density 
matrices instead of wave functions. 
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In classical Hamiltonian systems, the transition from regular motion to chaos 
is most clearly visible in the change of the phase-space structure: With increasing 
nonlinearity, regular tori successively dissolve in adjacent chaotic layers which 
grow in size and merge until the whole phase space is uniformly covered by a 
chaotic sea where the dynamics is locally hyperbolic and globally diffusive [3]. 
Research in quantum chaos has initially been concentrated on this limiting case 
of "hard chaos" , because the absence of structure in phase space facilitates the 
description. 

Closer to the generic situation, however, is the intermediate regime with an 
extremely intricate interweaving of regular and chaotic areas, as described by 
the Kolmogorov-Arnol'd-Moser (KAM) theorem, with self-similar hierarchies of 
regular islands. It is in this regime that we expect the most interesting, but at 
the same time least tractable, phenomena of chaos-coherence interplay to occur. 
A prominent example is chaotic tunneling, the coherent exchange of probability 
between symmetry-related regular regions that are separated dynamically by a 
chaotic layer, instead of a static potential barrier [4-20]. Chaotic tunneling comes 
about by the simultaneous action of classical nonlinear dynamics and quantum 
coherence. Tunneling is extremely sensitive to any disruption of coherence as 
it occurs due to the unavoidable coupling to the environment: In presence of 
dissipation, coherent tunneling becomes a transient that fades out on the way to 
an asymptotic state [21,22]. This is just one instance of the general rule that de- 
coherence tends to restore classical behaviour, other examples being the partial 
lifting, by dissipation, of the quantum suppression of chaos, and the appearance 
of quantum stationary states that show a close resemblance to corresponding 
classical strange attractors [2]. However, particularly for weak dissipation, more 
complicated cross effects occur, such as the strong modification of the decoher- 
ence time by chaotic tunneling. 

In this contribution, we investigate the mutual influence of chaotic tunneling 
and dissipation for a specific, but nevertheless generic case: a periodically forced 
bistable system. The quartic double well with a harmonic driving will serve as 
our working model. In Section 2 we introduce its Hamiltonian and the underlying 
symmetries. To provide the necessary background, we also briefiy review other 
important features of this system, in particular driven tunneling and its coherent 
suppression and modification in the presence of classical chaos without damping. 

Dealing with a driven system, its quantum dynamics is adequately analyzed 
in terms of the Floquet or quasienergy spectrum, also introduced in Section 2. 
The quasispectrum associated with chaotic tunneling exhibits a characteristic 
feature: Quasienergies of chaotic singlets frequently intersect tunnel doublets 
which are supported by regular tori. As the basic mechanism of chaotic tunneling 
we study, in Sections 3 and 4, the coherent and dissipative quantum dynamics 
in the vicinity of such singlet-doublet crossings. While in the coherent case the 
dynamics is well described in a three-state approximation, the coupling to the 
environment indirectly couples the three states to all other states. On the basis 
of numerical results for the full driven double well with dissipation, we reveal 
the limitations of the three-level approximation and identify additional features 
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of the full dynamics not covered by it. In particular, we consider the long-time 
asyniptotics and the phase-space structure associated with it. 

Also on the classical level, the presence of friction has profound consequences 
for the phase-space structure: Due to the net contraction of phase-space volume, 
stationary states are restricted to manifolds of lower dimensionality than the un- 
derlying phase space. Depending on friction strength and details of the system, 
this attractor may be consist of fixed points, of limit cycles, or, if the classical 
dynamics is chaotic, of a strange attractor with self-similar, fractal geometry. On 
a quantum level, the structures associated with classical attractors are smeared 
out on a scale h, yet leave clear traces in the asymptotic state of the corre- 
sponding dissipative quantum dynamics [23]. We study the classical-quantum 
correspondence of the asymptotic state in Section 5. 

2 The model 

We consider the quartic double well with a spatially homogeneous driving force, 
harmonic in time. It is defined by the Hamiltonian 

H{t)^HY)w+HFit), (1) 

Hpit) = Sxcos{f2t). (3) 

The potential term of the static bistable Hamiltonian i^DW possesses two minima 
at a; = ±xo, xo = (8-Eb/wWq)^/^, separated by a barrier of height Eb (cf. Fig. 1). 
The parameter loq denotes the (angular) frequency of small oscillations near the 
bottom of each well. Apart from mere scaling, the classical phase space of ^dw 
only depends on the presence or absence, and the signs, of the x"^ and the cc^ 
term. Besides that, it has no free parameter. This is obvious from the scaled 
form of the classical equations of motion, 

x=p, (4) 

p-=^x~^x^~Fcos{nt), (5) 

where the dimensionless quantities x, p and t are given by x/xq, p/muJoXo and 
uiot, respectively. The influence of the driving on the classical phase-space struc- 
ture is fully characterized by the rescaled amplitude and frequency of the driving, 

F = ^ ^ ^ n^ — . (6) 

y'SmujQEB ^0 

This implies that the classical dynamics is independent of the barrier height Eb- 
In the quantum-mechanical case, this is no longer true: The finite size of 
Planck's constant results in a finite number 
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Fig. 1. Sketch of the driven dou- 
ble well potential described by the 
time-dependent Hamiltonian (1) at 
two different phases. 
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of doublets with energy below the barrier top. The classical limit amounts to 
D — > CX3. This is evident from the classical scales for position, xq, and momen- 
tum, mujoxo, introduced above: The corresponding action scale is mcooXQ and 
therefore, the position-momentum uncertainty relation in the scaled phase space 
(x,p) reads 

AxAp>^, (8) 



where 



1 
8D 



(9) 



mutoXQ 

denotes the effective quantum of action. 

In the following, we restrict the driving amplitude to moderate values, such 
that the difference between the potential minima remains much smaller than the 
barrier height. This implies that the bistable character of the potential is never 
lost. 



2.1 Symmetries 

Time periodicity. The Hamiltonian (1) is P-pcriodic, with P = 27r/i7. As 
a consequence of this discrete time-translational invariance of H{x,p;t), the 
relevant generator of the quantum dynamics is the Floquet operator [24-28] 



U ^Tcxp -- / dtH{t) 



(10) 



where T denotes time ordering. According to the Floquet theorem, the adiabatic 
states of the system are the eigenstates of U. They can be written in the form 



\Mt)) 



''^""IMt)), 



(11) 
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with 

\c^^{t + P))^\c^^{t)). 
Expanded in these Floquet states, the propagator of the driven system reads 

U{t',t) = ^C--°(*'-*)/^|0„(t'))('^aW|- (12) 

a 

The associated eigenphases Cq,, referred to as quasienergies, come in classes, 
^a,n = Eq + nhf2, n ~ 0, ±1, ±2, . . .. This is suggested by a Fourier expansion of 
the \Mt)), 



\(f>a{t)) = y^|Ca,n)e 



%nt 



|ca,„)- p / dt|(/)„(i))e"^*. (13) 

^ Jo 

The index n counts the number of quanta in the driving field. Otherwise, the 
members of a class a are physically equivalent. Therefore, the quasienergy spec- 
trum can be reduced to a single "Brillouin zone", —fifi/2 < e < fifi/2. 

Since the quasienergies have the character of phases, they can be ordered 
only locally, not globally. A quantity that is defined on the full real axis and 
therefore does allow for a complete ordering, is the mean energy [22,27,28] 

1 /-^ 

Ec.= pJ dt(^„(i)|ff(t)|V'a(i)) = ((</>a(t)|//(i)|0a(i))). (14) 

It is related to the corresponding quasienergy by 

£^«=6a + ((<^«(i)|i?i^|0a(t))), (15) 

where the outer angle brackets denote the time average over one period of the 
driving, as indicated by Eq. (14). The second term on the right-hand side plays 
the role of a geometric phase accumulated over this period [22,29]. Without the 
driving, E^ — £«, as it should be. By inserting the Fourier expansion (13), the 
mean energy takes the form 

Ea ^^i'^a+ nflQ) {Ca,n\Ca,n)- (16) 

n 

It shows that the nth Floquet channel gives a contribution ea + nTiil to the mean 
energy, weighted by the Fourier coefficient {Ca^n\ca,n) [28]. 

Quasienergies and Floquet states are obtained numerically by solving the 
matrix eigenvalue equation [24, 27, 28] 

2_^/_^'^n,k;n',k'Cn',k'—^Cn,ki (17) 
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equivalent to the time-dependent Schrodinger equation. It is derived by inserting 
the eigenstates (11) into the Schrodinger equation, Fourier expanding, and using 
the representation in the eigenbasis of the unperturbed Hamiltonian, Ho\Wk) = 
EkWk)- We introduced the abbreviations 

T-Ln 



•■,n',k' 


^{Ek^ 


nfif2)5n- 


-n'Sk- 


-k' 


+ 


1 

2 


SxkM' 


{5„- 


-1- 


-n' 


+ Sn+1- 


-«'), 


(18) 


Cn,k 


= {^k\Cr. 
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(20) 



Time-reversal symmetry. The energy eigenfunctions of an autonomous Ham- 
iltonian with time-reversal symmetry, 

T : X ^ X, p ^ —p, t ^ —t (21) 

can be chosen as real [30,31]. Time-reversal invariance is generally broken by 
a magnetic field or by an explicit time-dependence of the Hamiltonian. How- 
ever, for the sinusoidal shape of the driving together with the initial phase 
chosen above, T invariance is retained and the Schrodinger operator 7i(i) = 
H{t) — \Tidt obeys 7i(t) = Ti*{—t). If now (j>{x,t) is a Floquet state in position 
representation with quasienergy e, then 0*(x, —t) also is a Floquet state with the 
same quasienergy. This means that we can always find a linear combination of 
these Floquet states such that (j>(x, t) = 4>*{x, —t), or in the frequency domain, 
(j){x,n) — (l)*(x,n), i.e., the Fourier coefficients of the Floquet states can be 
chosen real. 



Generalized parity. The invariance of ifow under parity P: x ^ — x, p -^ —p, 
t —* t is destroyed by any spatially constant driving force. With the above choice 
of Hp{t), however, a more general, dynamical symmetry remains [32-34]. It is 
defined by the operation 

Pp : X ^ -X, p^ -p, t^t + P/2 (22) 

and represents a generalized parity acting in the extended phase space spanned 
by X, p, and phase, i.e., time tniodP. While such a discrete symmetry is of 
minor importance in classical physics, its influence on the quantum mechanical 
quasispectrum {ea{F)} is profound: It devides the Hilbert space in an even 
and an odd sector, thus allowing for a classification of the Floquet states as 
even or odd. Quasiencrgies from different symmetry classes may intersect, while 
quasienergies with the same symmetry typically form avoided crossings [31]. 
The fact that Pp acts in the phase space extended by time imodP, results in a 
particularity: If, e.g., |(/)(t)) is an even Floquet state, then exp(iJ7t)|(/)(t)) is odd, 
and vive versa. Thus, two equivalent Floquet states from neighboring Brillouin 
zones have opposite generalized parity. This means that a classification of the 
corresponding solutions of the Schrodinger equation, Itpit)) = exp{—[et/h)\4>{t)), 
as even or odd is meaningful only with respect to a given Brillouin zone. 
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The invariancc of the system under Pp is also of considerable help in the 
numerical treatment of the Floquet matrix (18) [11, 12]. To obtain a complete 
set of Floquet states, it is sufficient to compute all eigenvectors of the Floquet 
Hamiltonian in the even subspace whose eigenvalues lie in the first two Brillouin 
zones. The even Floquet states are given by the eigenvectors of He from the 
first Brillouin zone; the odd Floquet states are obtained by shifting the (even) 
ones from the second to the first Brillouin zone, which changes their generalized 
parity. Thus, in the even subspace, we have to diagonalize the matrix 



/. 



He 
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■ ■ Ee + 2hn 
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(23) 



For the same number of Floquet channels, it has only half the dimension of the 
original Floquet matrix (18). The entries in He are themselves blocks of infinite 
dimension, in principle. They read explicitly 



Ee = 



/Eo ■■■\ 

E2 ■■ 

£^4 • • 



V 



Eo = 



(El •••\ 

£^3 ••• 
£5 • • • 



•7 
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(24) 
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7 



Here, -Ec, E^, represent the undriven Hamiltonian, and Xoo, -^^oo the driving 
field Hi = Sx/2, decomposed in the even and odd eigenstates \'Fk) of H^w, 
with Ek denoting an eigenvalue of TJdWi ^-^d Xk^k' a matrix element of the 
position operator, see Eq. (20). 



2.2 Tunneling, driving, and dissipation 

With the driving Hp (t) switched off, the classical phase space generated by -ffDW 
exhibits the constituent features of a bistable Hamiltonian system. A separatrix 
ai E = forms the border between two sets of trajectories: One set, with E < 0, 
comes in symmetry-related pairs, each partner of which oscillates in either one of 
the two potential minima. The other set consists of unpaired, spatially symmetric 
trajectories, with E > 0, which encircle both wells. 
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Torus quantization of the integrable undrivcn double well, Eq. (2), implies 
a simple qualitative picture of its eigenstates: The unpaired tori correspond 
to singlets with positive energy, whereas the symmetry-related pairs below the 
top of the barrier correspond to degenerate pairs of eigenstates. Due to the 
almost harmonic shape of the potential near its minima, neighboring pairs are 
separated in energy approximately by fiujQ. Exact quantization, however, predicts 
that the partners of these pairs have small but finite overlap. Therefore, the true 
eigenstates come in doublets, each of which consists of an even and an odd state, 
\^n) and |^~), respectively. The energies of the nth doublet are separated by a 
finite tunnel splitting Z\„. We can always choose the global relative phase such 
that the superpositions 

|#^>L) = -^ (!<?+) ± \<P-)) (26) 

are localized in the right and the left well, respectively. As time evolves, the 
states |<?^), |<?~) acquire a relative phase exp(— iZ\„i/?i) and |^^), |<?^) are 
transformed into one another after a time nh/An- Thus, the particle tunnels 
forth and back between the wells with a frequency An/h. This introduces an 
additional, purely quantum-mechanical frequency scale, the tunnel rate Aq/H 
of a particle residing in the ground-state doublet. Typically, tunnel rates are 
extremely small compared to the frequencies of the classical dynamics, all the 
more in the semiclassical regime we are interested in. 

A driving of the form (3), even if its influence on the classical phase space 
is minor, can entail significant consequences for tunneling: It may enlarge the 
tunnel rate by orders of magnitude or even suppress tunneling altogether. For 
adiabatically slow driving, i? ^ Ao/h, tunneling is governed by the instanta- 
neous tunnel splitting, which is always larger than its unperturbed value Aq 
and results in an enhancement of the tunneling rate [33]. If the driving is faster, 
Ao/h ^ /2 ^ Wo, cf. Fig. 2, the opposite holds true: The relevant time scale is 
now given by the inverse of the quasienergy splitting of the ground-state doublet 
h/\ei — e()\. It has been found [33,35] that in this case, for finite driving amplitude, 
I El ~ Col < ^0- Thus tunneling is always decelerated. Where the quasienergies 
of the ground-state doublet (which are of different generalized parity) intersect 
as a function of F, the splitting vanishes and tunneling is brought to a complete 
standstill by the purely coherent influence of the driving [32] . 

The small energy scales associated with tunneling make it extremely sensi- 
tive to any loss of coherence. As a consequence, the symmetry underlying the 
formation of tunnel doublets is generally broken, and an additional energy scale 
is introduced, the effective finite width attained by each discrete level. Tunneling 
and related coherence phenomena thus fade out on a time scale tdocoh- In gen- 
eral, this time scale gets shorter for higher temperatures, reflecting the growth 
of the transition rates (53) [36]. However, there exist counterintuitive effects: For 
driven tunneling in the vicinity of an exact crossing of the ground-state doublet, 
the coherent suppression of tunneling [22, 32, 33] can be stabilized with higher 
temperatures [37-39] until levels outside the doublet start to play a role. 
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Fig. 2. Tunneling plienomena and tlie according appropriate levels of description for the 
non-dissipative driven double- well potential, Eq. (1). The bars depict the corresponding 
regimes of the driving frequency Q. See Section 2 for a detailed discussion. 



So far, we have considered only driving frequencies much smaller than the 
frequency scale wq of the relevant classical resonances, i.e., a parameter regime 
where classical motion is predominantly regular. In this regime, coherent tunnel- 
ing is well described within a two-state approximation [33,35]. In the dissipative 
case, however, a two-state approximation must fail for temperatures k^T ^ Tiujq, 
where thermal activation to higher doublets becomes relevant. 

2.3 The onset of chaos 

Driving with a frequency i7 « wq affects also the dynamics of the classical 
bistable system, as small oscillations near the bottom of the wells become reso- 
nant with the driving and classical chaos comes into play (cf. Fig. 2). In a quan- 
tum description, this amounts to resonant multiple excitation of inter-doublet 
transitions until levels near the top of the barrier are significantly populated. 

In this frequency regime, switching on the driving has two principal conse- 
quences for the classical dynamics: The separatrix is destroyed as a closed curve 
and replaced by a homoclinic tangle [40] of stable and unstable manifolds. A 
chaotic layer forms in the vicinity and with the topology of the former separa- 
trix (cf. Fig. 6, below). This opens the way for diffusive transport between the 
two potential wells. Due to the nonlincarity of the potential, there is an infinite 
set of resonances of the driving with the unperturbed motion, both inside as well 
as outside the wells [41,42]. Since the period of the unperturbed, closed trajecto- 
ries diverges for i? ^ 0, both from below and above, the resonances accumulate 
towards the separatrix of the unperturbed system. By its large phase-space area, 
the first resonance (the one for which the periods of the driving and of the unper- 
turbed oscillation are in a 1-to-l ratio) is prominent among the others and soon 
(in terms of increasing amplitude F) exceeds in size the "order-zero" regular 
areas near the bottom of each well [11]. 
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Both major tendencies in the evolution of the classical phase space, extension 
of the chaotic layer and growth of the first resonance, leave their specific traces in 
the quasienergy spectrum. The tunnel doublets characterizing the unperturbed 
spectrum for _E < pertain to states located on pairs of symmetry-related 
quantizing tori in the regular regions within the wells. With increasing size of 
the chaotic layer, the quantizing tori one by one resolve in the chaotic sea. 
The corresponding doublets disappear as distinct structures in the spectrum 
as they attain a splitting of the same order as the mean level separation. The 
gradual widening of the doublets proceeds as a smooth function of the driving 
amplitude [11, 12], which roughly obeys a power law [43-45]. As soon as a pair 
of states is no longer supported by any torus-like manifold, including fractal [46] 
and vague tori [47], the corresponding eigenvalues detach themselves from the 
regular ladder to which they formerly belonged. They can then fluctuate freely 
in the spectrum and thereby "collide" with other chaotic singlets or regular 
doublets. 

The appearance of a regular region, large enough to accommodate several 
eigenstates, around the first resonance introduces a second ladder of doublets 
into the spectrum. Size and shape of the first resonance vary in a way different 
from the main regular region. The corresponding doublet ladder therefore moves 
in the spectrum independently of the doublets that pertain to the main regular 
region, and of the chaotic singlets. This gives rise to additional singlet-doublet 
and even to doublet-doublet encounters. 



3 Chaotic tunneling near singlet-doublet crossings 



Near a crossing, level separations deviate vastly, in both directions, from the 
typical tunnel splitting (cf. Fig. 8, below). This is reflected in time-domain phe- 
nomena ranging from the suppression of tunneling to a strong increase in its 
rate and to comphcated quantum beats [13-15]. Singlet-doublet crossings, in 
turn, drastically change the non-dissipative quasienergy scales and replace the 
two-level by a three-level structure. As a consequence, the familiar way tunnel- 
ing fades out in the presence of dissipation is also signiflcantly altered. Near a 
crossing, the coherent dynamics can last much longer than for the unperturbed 
doublet, despite the presence of the same dissipation as outside the crossing, 
establishing "chaos-induced coherence." Depending on temperature, it can also 
be destroyed on a much shorter time scale. 

For the parameters chosen in our numerical studies, higher resonances are 
negligible in size. Accordingly, the "coastal strip" between the chaotic layer along 
the former separatrix and the regular regions within and outside the wells, formed 
by hierarchies of regular islands around higher resonances, remains narrow (cf. 
Fig. 6). For the tunneling dynamics, the role of states located in the border 
region [17, 18] is therefore not significant in our case. 
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Fig. 3. Possible configurations of qua- 
sienergy crossings between a cliaotic 
singlet and a regular doublet. Differ- 
ent line types signify different parity. 
See Section 3.1 for the labeling of the 
levels. Note that only for configura- 
tions (a),(b), the order of the regular 
doublet is restored in passing through 
the crossing. In configurations (c),(d), 
it is reversed. 



3.1 Three- level crossings 

Among the various types of quasienergy crossings that occur according to the 
above scenario, those involving a regular doublet and a chaotic singlet are the 
most common. In order to give a quantitative account of such crossings and 
the associated coherent dynamics, and for later reference in the context of the 
incoherent dynamics, we shall now discuss them in terms of a simple three-state 
model, devised much in the spirit of Ref . [7] . 

Far to the left of the crossing, we expect the following situation: There is a 
doublet of Floquet states 



|V+(t))=c 



= „-ie+i/fi 



\^nt)), 



lV'r(i)) = c 



= e-'K++^)*Au- 



\<f'7{t)). 



(27) 
(28) 



with even (superscript +) and odd (— ) generalized parity, respectively, resid- 
ing on a pair of quantizing tori in one of the regular (subscript r) regions. We 
have assumed the quasienergy splitting Z\ = e~ — e+ (as opposed to the unper- 
turbed splitting) to be positive. The global relative phase is chosen such that 
the superpositions 



l0R,L(t)) 



1 

V2 



(i0+(t))±i0r(i))) 



(29) 



are localized in the right and the left well, respectively, and tunnel back and 
forth with a frequency A/h. 

As the third player, we introduce a Floquet state 



IC(i)>=c 



= o-i(£ ++-4+'4c)i/fi| A- 



IC(t)>, 



(30) 



located mainly in the chaotic (subscript c) layer, so that its time-periodic part 
|(/)~(t)) contains a large number of harmonics. Without loss of generality, its 
parity is fixed to be odd. For the quasienergy, we assume that e^T = e++Zi+Z\c = 
e~ + Ac, where |Z\c| can be regarded as a measure of the distance from the 
crossing. 
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The structure of the classical phase space then implies that the mean energy 
of the chaotic state should be close to the top of the barrier and far above that 
of the doublet. We assume, like for the quasienergics, a small splitting of the 
mean energies pertaining to the regular doublet, \E:^ — _E+| ^ E^ — Ef. 

In order to model an avoided crossing between \(t>^) and Ii/i^t), we suppose 
that there is a non- vanishing fixed matrix clement 

b^{{cl,;\Huw\<Pc))>0. (31) 

For the singlet-doublet crossings under study, we typically find that Z\ <C 6 <C 
hi7. Neglecting the coupling with all other states, we model the system by the 
three-state (subscript 3s) Floquet Hamiltonian 

/O \ 
Has - e+ + Z\ b (32) 

\0 b A + aJ 

in the three-dimensional Hilbert space spanned by {|0+(i)), |07(t)), \4>^{t))}. Its 
Floquet states are 

|^+(t)) = |0+(t)), 

l^iit)) - (|0r(i)>cos/3- |C(i))sin/3) , (33) 

I02"(i)> - (|0r(i)>sin/3+lC(t))cos/3). 

with quasienergics 

1 1 

2 '='^2 



6+-6+, e^^,^et+A + -A,T7;V^+^, (34) 



and mean energies, neglecting contributions of the matrix element b, 

(35) 



1+ 

'0 


= e;, 


'1 


= E- cos^ 13 + E- sin^ /?, 


'2 


= E- sin^ 13 + E- cos^ /3. 



The angle (3 describes the mixing between the Floquet states |(^~) and \4>~) 
and is an alternative measure of the distance to the avoided crossing. By diago- 
nalizing the Hamiltonian (32), we obtain 

2/3 = arctan (-r-j, Q < [3 <^. (36) 

For (3 -^ 7r/2, corresponding to —A^ ^ 6, we retain the situation far left of the 
crossing, as outlined above, with \cf)'^) w \(j)c)^ l'/'2^) ~ \'t'7)- To the far right of 
the crossing, i.e., for /3 ^- or Ac » b, the exact eigenstates \(j)i) and |(/)^) have 
interchanged their phase-space structure [13-15]. Here, we have \(t)i) ~ \4>7) 
and 10^) « \4'c)- The mean energy is essentially determined by this phase-space 
structure, so that there is also an exchange of E^ and E2 in an exact crossing, 
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Fig. 4. A singlet-doublet crossing, according to the three-state model (32), in terms of 
the quasienergies (a) and the mean energies (b) as functions of the coupling parameter 
Z\c/&. Energies for a corresponding exact crossing (i.e., with the crossing states uncou- 
pled) are marked by dotted lines, the energies in the presence of coupling by full and 
dashed lines for even and odd states, respectively. 
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Fig. 5. Singlet-doublet crossing found numerically for the driven double well, Eq. (1), 
at _D = 4 and i? = 0.982 a;o, in terms of the dependence of the quasienergies (a) and the 
mean energies (b) on the driving amplitude F . Values of the driving amplitude used 
in Fig. 9 are marked by dotted vertical lines. Full and dashed lines indicate energies 
of even and odd states, respectively. Bold lines give the mean energies of the chaotic 
singlet and the ground-state doublet depicted in panel (a). 



14 



Peter Hanggi, Sigmund Kohler, and Thomas Dittrich 



0.5 



3 0.0 



-0.5 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

x/xq 

Fig. 6. Classical stroboscopic phase space portrait, aX t = 27rn/J7, of the harmonically 
driven quartic double well, Eq. (1). The driving parameters F — 0.015, O — 0.982 ojo, 
are at the center of the singlet-doublet crossing shown in Fig. 5. 



0.5 



3 0.0 



-0.5 



(a) 





0.5 



3 0.0 



-0.5 




-1.5 -1.0 -0.5 0.0 0.5 1.0 
x/xo 



1.5 



-1.5 -1.0 -0.5 0.0 0.5 1.0 
x/xo 



1.5 



Fig. 7. Contour plots of the Husimi functions for the Floquet states \(l>i) ~ \4'7) 
(a) and |(^2^) ~ I^^T) (b) of the harmonically driven quartic double well, Eq. (1), at 
stroboscopic times t — nP. The driving parameters F = 0.014, i? = 0.982 tJo, are in 
sufficient distance to the singlet-doublet crossing such that the admixture from the 
chaotic singlet state is negligible. The rectangle in the lower left corner has the size of 
the effective quantum of action Ticft- 



cf. Eq. (35), while Eq remains unaffected (Fig. 4b). The quasienergies e^ and ej~ 
must intersect close to the avoided crossing of e^ and e^ (Fig. 4a). Their crossing 
is exact, since they pertain to states with opposite parity (cf. Fig. 3a, b). 

In order to illustrate the above three-state model and to demonstrate its 
adequacy, we have numerically studied a singlet-doublet crossing that occurs 
for the double- well potential, Eq. (1), with 13 = 4, at a driving frequency J? = 
0.982 ciJo and amplitude F — 0.015029 (Fig. 5). The phase-space structure of the 
participating Floquet states (Figs. 6, 7) meets the assumptions of our three-state 
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Fig. 8. Splitting of the low- 
est doublets for D — 4 and 
n = 0.982 Wo- The arrows indi- 
cate the locations of the exact 
and the avoided crossing within 
a three- level crossing of the type 
sketched in Fig. 3a. 



theory. A comparison of the appropriately scaled three-state theory (Fig. 4) with 
this real singlet-doublet crossing (Fig. 5) shows satisfactory agreement. Note 
that in the real crossing, the quasienergy of the chaotic singlet decreases as 
a function of F, so that the exact crossing occurs to the left of the avoided 
one. This numerical example also shows a deficiency of the idealized three-state 
model. Following the global tendency of widening of the splittings with increasing 
driving amplitude [11,44,45], it may happen that even far away from a crossing, 
the doublet sphtting does not return to its value on the opposite side (see Fig. 8). 
It is even possible that an exact crossing of ej and ej" never takes place in the 
vicinity of the crossing. In that case, the relation of the quasienergies in the 
doublet gets reversed via the crossing (Fig. 3c, d). Nevertheless, the three-state 
scenario captures the essential features. 

To study the dynamics of the tunneling process, we focus on the state 



m)) = -^ (c- 



[t/h 



I4W> 



-'/^\^- 



{t))cosP + e' 



r*/^|02-(t))sin/3 



(37) 

It is constructed such that at i = 0, it corresponds to the decomposition of J^r) 
in the basis (33) at finite distance from the crossing. Therefore, it is initially 
localized in the regular region in the right well and follows the time evolution 
under the Hamiltonian (32). From Eqs. (29), (33), we find the probabilities for 
its evolving into |0r), |0l), or \(l)c), respectively, to be 



PnMt) 




{t)mm' 



i^i 



4)t 



cos P + cos 



(^2 



COS 



, (ei - €2 )t -I 

5 : — i 



p,{t) = \{Mt)mtw = 



1 — COS 




-^sin2/3 



(38) 



cos^/3sin /3. 
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We discuss the coherent dynamics of the three-state model for different dis- 
tances to the crossing and illustrate it by numerical results for the real crossing 
introduced above. 

At sufficient distance from the crossing, there is only little mixing between the 
regular and the chaotic states, i.e., sin/3 ^ 1 or cos/3 ^ 1. The tunneling process 
then follows the familiar two-state dynamics involving only \4>^) and {(j)^), with 
tunnel frequency A/h (Fig. 9a). Close to the avoided crossing, cos/3 and sin/3 
are of the same order of magnitude, and \4>i), \4>2) become very similar to one 
another. Both now have support in the chaotic layer as well as in the symmetry- 
related regular regions, they are of a hybrid nature. Here, the tunneling involves 
all the three states and must be described at least by a three-level system. The 
exchange of probability between the two regular regions proceeds via a "stop- 
over" in the chaotic region [7,8,13-15]. 

The three quasienergy differences that determine the time scales of this pro- 
cess are in general all different, leading to complicated beats (Fig. 9b). However, 
for Ac = — 2Z\, the two quasienergies e^ — e'^ and e^ — €2 are degenerate. At 
this point, the center of the crossing, the number of different frequencies in the 
three-level dynamics reduces to two again. This restores the familiar coherent 
tunneling in the sense that there is again a simple periodic exchange of prob- 
ability between the regular regions [13-15]. However, the rate is much larger 
if compared to the situation far off the crossing, and the chaotic region is now 
temporarily populated during each probability transfer, twice per tunneling cycle 
(Fig. 9c). 

4 Incoherent quantum dynamics 

4.1 Master equation 

System-bath model. To achieve a microscopic model of dissipation, we couple 
the system ( 1 ) bilinear ly to a bath of non- interacting harmonic oscillators [48 , 49] . 
The total Hamiltonian of system and bath is then given by 

H{t) = HuMt) + E f ^ + ^-l U - ^-)1 • (39) 

The position x of the system is coupled, with coupling strength gj^, to an ensemble 
of oscillators with masses TOj^, frequencies w^, momenta p^, and coordinates x^. 
The bath is fully characterized by the spectral density of the coupling energy, 

00 2 

J(c.)=^^-^^^(^-c.,). (40) 

u—l 

For the time evolution we choose an initial condition of the Fcynman- Vernon 
type: at t — to, the bath is in thermal equilibrium and uncorrelated to the 
system, i.e., 

P{to) =^ Psih) <» PB,cq, (41) 



(i'i^'>'"i;v^ 
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Fig. 9. Stroboscopic time evolution 
of a state initially localized in the 
right well, in the vicinity of the 
singlet-doublet crossing shown in 
Fig. 5, in terms of the probabili- 
ties to be in the right well (which 
here is identical to the return prob- 
ability, marked by full lines), in 
the reflected state in the left well 
(dashed), or in the chaotic state 
\tp^} (dotted). Parameter values are 
as in Fig. 5, and F = 0.0145 (a), 
0.0149 (b), 0.015029 (c). 



where pB.eq = cxp(— /3i/B)/trB cxp(— /^i/e) is the canonical ensemble of the bath 
and 1//3 = kBT. 

Due to the bilinearity of the system-bath coupling, one can always eliminate 
the bath variables to get an exact, closed integro-differential equation for the 
reduced density matrix ps(t) = trBp(t). It describes the dynamics of the central 
system, subject to dissipation [50]. 



Born-Markov approximation. In most cases, however, the integro-differen- 
tial equation for ps{t) can be solved only approximately. In particular, in the 
limit of weak coupling. 



7 < kBT/h, 

7< ka - ea'\/h, 



(42) 
(43) 



it is possible to truncate the time-dependent perturbation expansion in the 
system-bath interaction after the second-order term. The quantity 7, to be 
defined below, denotes the effective damping of the dissipative system, and 
l^a ~ Ca'l/^ ^^^ t^^ transition frequencies of the central system. In the present 
case, the central system is understood to include the driving [51-54], so that 
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the transition frequencies are given by quasienergy differences. The autocorre- 
lations of the bath decay on a time scale Ti/k^T and thus in the present limit, 
instantaneously on the time scale I/7 of the system correlations. With the initial 
preparation (41), the equation of motion for the reduced density matrix in this 
approximation is given by [54] 

Ps{t) - -i [Hs{t),ps{t)] + 4r / ^^J{^)ntb{^Lo) 

/•OO 

x/ dT{e''^^[£{t-T,t)ps{t),x]+R.c.), (44) 

Jo 

where x{t' , t) denotes the position operator in the interaction picture defined by 

x{t',t)^U\t',t)xU{t',t), (45) 

with U{t',t), the propagator of the conservative driven double well, given in 
Eq. (12). 'H.c' means 'Hermitian conjugate', and 

"th(e) = g,/fc3T _ 1 = -"th(-e) - 1 (46) 

is the thermal occupation of the bath oscillator with energy e. To achieve a more 
compact notation, we require J{—uj) = —J{uj). In the following, we shall restrict 
ourselves to an Ohmic bath, J{lu) = Twyuj. This defines the effective damping 
constant 7. 

We use the time-periodic components \4>o,{t)) of the Floquet states as a basis 
to expand the density operator, Eq. (44). Expressing the matrix elements 

Xcpit) = {ct>^{t)\x\4>p{t)) (47) 

of the position operator by their Fourier coefficients 

^a/3(t)=^e"^^*X„^,„, (49) 

n 

yields the equation of motion for the elements p^p of the reduced density matrix 
PS [37,39,52,54], 

Po.p{t) = ^^{Mt)\psit)\Mt)) 

= "■^(^" -ef3)pai3{t) 

+ 2_^ [Naa'.nXaa' ,nPa'f3'Xf3'0,n' 
a' P'nn' 

- 7V„.;5.,„X„,.,„.X,„.,„p/3./3)e'("+"')'"* + H.c. . (50) 
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The coefficients of this differential equation are periodic in time with the period 
of the driving. The Nap^n are given by 

iVa/3,n = iV(e„ -ep + nhn), N{e) = ^nth(e). (51) 

For e 3> kBT, N{e) approaches zero. 

Since the position operator x is odd under Pp (cf. Eq. (22)), the master 
equations (44) and (50) are invariant under Pp. Therefore, both the conservative 
and the dissipative dynamics preserve the parity of the operator |(/)q,)(0^|. If 
|(/)q,) and \(j)fj) belong to the same parity class, it is even, and odd otherwise. In 
particular, the projectors \(j)a){4'a\ and thus all density matrices diagonal in the 
Floquet basis are even under Pp. 



Rotating-wave approximation. Assuming that dissipative effects are rele- 
vant only on a time scale much larger than the period P of the driving, we 
average the coefficients of the master equation (50) over P to obtain the equa- 
tion of motion 

Paf3{t) = ^^(ea ~ f-l3)Pap{t) + ^ Cafl^a' f3' Pa' f3' , (52) 

^ a'l3' 

with the time-independent dissipative part 

Cap,a'l3' — 2_^ {Naa'.n + ^/3/3',n) Xaa' ,nX fj' fj^-n 
n 

f3",n 

— Oaa' / ^ Naiip' ,nXpia" ,-nXa"l3,n- (53) 

a" n 

This step amounts to a rotating-wave approximation which is, however, less re- 
strictive than the one introduced in [51,52] where dissipative effects are averaged 
over the generally longer time scale maxa^^_„(27r?i/(eQ — ^p + nfiQ)). 

4.2 Chaos-assisted dissipative tunneling 

The crucial effect of dissipation on a quantum system is the disruption of coher- 
ence: a coherent superposition evolves into an incoherent mixture. Thus, phe- 
nomena based on coherence, such as tunneling, are rendered transients that fade 
out on a finite time scale idccoh- In general, for driven tunneling in the weakly 
damped regime, this time scale gets shorter for higher temperatures, as transition 
rates grow [36]. However, in the vicinity of an exact crossing of the ground-state 
quasienergies, the coherent suppression of tunneling [22,32,33] can be stabilized 
with higher temperatures [37-39] and increasing friction [55,56] until levels out- 
side the doublet start to play a role. We have studied dissipative chaos-assisted 
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Fig. 10. Occupation probabilities as in Fig. 9a,c, but in the presence of dissipation. The 
dash-dotted hne shows the time evolution of trp^. The parameter values are D = 4, 
f2 = 0.982 Wo, 7 = W'^uJo, fceT = IQ-^ftt^o, and F = 0.0145 (a), 0.015029 (b). The 
inset in (a) is a blow up of the rectangle in the upper left corner of that panel. 
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Fig. 11. Time evolution of the return 
probability Pr (full line) and the co- 
herence function trp^ (dash-dotted) 
during loss and regain of coher- 
ence. The parameter values are as in 
Fig. 10b. 



tunneling, at the particular real singlet-doublet crossing introduced in Sec. 3.1 
(see Fig. 5). The time evolution has been computed numerically by iterating 
the dissipative quantum map for the improved master equation in moderate 
rotating- wave approximation, Eq. (52). As initial condition, we have chosen the 
density operator p{{]) = |(/'r)(0r|, a pure state located in the right well. 

In the vicinity of a singlet-doublet crossing, the tunnel splitting increases 
significantly — the essence of chaos- assisted tunneling. During the tunneling, the 
chaotic singlet \(j)c) becomes populated periodically with frequency \€2^e.'^\/h, cf. 
Eq. (38) and Fig. 9. The high mean energy of this singlet results in an enhanced 
decay of coherence at times when it is well populated (Fig. 10). For the relaxation 
towards the asymptotic state, also the slower transitions within doublets are 
relevant. Therefore, the corresponding time scale troiax can be much larger than 
idocoh (Fig. 11). 

To obtain quantitative estimates for the dissipative time scales, we approxi- 
mate tdecoh by the decay rate of trp^, as a measure of coherence, averaged over 
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Fig. 12. Time scales of the decay of the coherence measure trp^ (a) and of the relax- 
ation towards the asymptotic solution (b) near the singlet-doublet crossing. Near the 
exact crossing (F « 0.013, full vertical line) coherence is stabilized, whereas at the cen- 
ter of the avoided crossing {F ^ 0.015, dashed vertical line) the decay of coherence is 
accelerated. The parameter values are D — 4, f2 — 0.982 loq, 7 = 10~^uJo, temperature 
as given in the legend. 



a time t 



P' 



^dccoh 



dt'-^trp'if) 
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"-p ^0 
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Because of the stepwise loss of coherence (Fig. 10), we have chosen the propa- 
gation time tp as an rifold multiple of the duration 27r?i/|e^ ~ ^i\ "^f the chaotic 
beats. For this procedure to be meaningful, n should be so large that the coher- 
ence decays substantially during the time tp (in our numerical studies to a value 
of approximately 0.9). The time scale ireiax of the approach to the asymptotic 
state is given by the reciprocal of the smallest real part of the eigenvalues of the 
dissipative kernel. 

Outside the singlet-doublet crossing we find that the decay of coherence and 
the relaxation take place on roughly the same time scale (Fig. 12). At F ^ 0.013, 
the chaotic singlet induces an exact crossing of the ground-state quasienergies 
(see Fig. 8), resulting in a stabilization of coherence with increasing temperature. 
At the center of the avoided crossing, the decay of coherence becomes much faster 
and is essentially independent of temperature. This indicates that transitions 
from states with mean energy far above the ground state play a crucial role. 



4.3 Asymptotic state 

As the dynamics described by the master equation (52) is dissipative, it converges 
in the long-time limit to an asymptotic state Poo{t)- In general, this attractor 
remains time dependent but shares the symmetries of the central system, i.e. 
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here, periodicity and generalized parity. However, the coefhcients (53) of the 
master equation for the matrix elements Pap^ valid within a moderate rotating- 
wave approximation, are time independent and so the asymptotic solution also 
is. The explicit time dependence of the attractor has been effectively eliminated 
by representing it in the Floquet basis and introducing a mild rotating-wave 
approximation. 

To gain some qualitative insight into the asymptotic solution, we focus on 
the diagonal elements 



*^aQL.a' a' / , -^ *aa',n |^aQ',n I i (X f^ (X ^ V / 



of the dissipative kernel. They give the rates of direct transitions from \4>a') to 
|(/>Q,). Within a full rotating-wave approximation [51,52], these were the only non- 
vanishing contributions to the master equation to affect the diagonal elements 
Paa of the density matrix. 

In the case of zero driving amplitude, the Floquet states \(t>a) reduce to the 
eigenstates of the undriven Hamiltonian i^DW- The only non- vanishing Fourier 
component is then Icq^o), and the quasienergies e^ reduce to the corresponding 
eigenenergies Ea- Thus Laa,a'a' only consists of a single term proportional to 
N{ea — ea')- It describes two kinds of thermal transitions: decay to states with 
lower energy and, if the energy difference is less than fceT, thermal activation 
to states with higher energy. The ratio of the direct transitions forth and back 
then reads 

*^aa,a' a' I ^a ^a' \ /r-^\ 

= cxp T^T^ • (57) 



We have detailed balance and therefore the steady-state solution is 

p„„.((X3)^e-'°/'=^^(5„„.. (58) 

In particular, the occupation probability decays monotonically with the energy 
of the eigenstates. In the limit /cbT -^ 0, the system tends to occupy the ground 
state only. 

For a strong driving, each Floquet state |(/)q,) contains a large number of 
Fourier components and Caa.a'a' is given by a sum over contributions with 
quasienergies ea — ea' + nfiil. Thus decay to states with "higher" quasienergy 
(recall that quasienergies do not allow for a global ordering) becomes possible 
due to terms with n < 0. Physically, it amounts to an incoherent transition 
under absorption of driving-field quanta. Correspondingly, the system tends to 
occupy Floquet states comprising many Fourier components with low index n. 
According to Eq. (16), these states have a low mean energy. 

The effects under study are found for a driving with a frequency of the 
order of unity. Thus for a quasienergy doublet, not close to a crossing, we have 
\ea ~ Eq'I ^ ^^1 and Ca'a'.aa IS dominated by contributions with n < 0, where 
the splitting has no significant influence. However, up to the tunnel splitting, 
the two partners in the quasienergy doublet are almost identical. Therefore, with 
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Fig. 13. Occupation probability paa of the Floquet states |0a) in the long-time limit. 
The parameter values are D — A, f2 = 0.982 tJo, 7 = 10^^ci;o, and F — 0.013 (a), 
0.015029 (b), temperature as given in the legend. 



respect to dissipation, both should behave similarly. In particular, one expects an 
equal population of the doublets even in the limit of zero temperature (Fig. 13a), 
in contrast to the time-independent case. 

In the vicinity of a singlet-doublet crossing the situation is more subtle. Here, 
the odd partner, say, of the doublet mixes with a chaotic singlet, cf. Eq. (33), 
and thus acquires components with higher energy. Due to the high mean energy 
E~ of the chaotic singlet, close to the top of the barrier, the decay back to 
the ground state can also proceed indirectly via other states with mean energy 
below E~ . Thus \(pi) and |0^) are depleted and mainly |0q ) will be populated. 
However, if the temperature is significantly above the splitting 26 at the avoided 
crossing, thermal activation from |0q ) to \4'i2)j accompanied by depletion via 
the states below E^ , becomes possible. Thus asymptotically, all these states 
become populated in a steady closed flow (Fig. 13b). The long-time limit of the 
corresponding classical dynamics converges to one of two limit cycles, each of 
which is located close to one of the potential minima. In a stroboscopic map they 
correspond to two isolated fixed points. This behavior is qualitatively different 
from the asymptotic limit of the dissipative quantum dynamics near the center 
of the crossing and shows that the occupation of levels outside the singlet and 
the doublet for i — > cx3 is a pure quantum effect. 

Important global characteristics of the asymptotic state, measuring its de- 
gree of spreading over phase space, are the Shannon entropy S — — tr (poc Inpoo) 
or, alternatively, trp^. The latter gives approximately the inverse number of in- 
coherently occupied states and can be considered an "incoherent inverse partici- 
pation ratio" [57]. It equals unity only if the attraetor is a pure state. According 
to the above scenario, we expect trp^ to assume the value 1/2, in a regime with 
strong driving but preserved doublet structure, reflecting the incoherent popu- 
lation of the ground-state doublet. In the vicinity of the singlet-doublet crossing 
where the doublet structure is dissolved, its value should be close to unity for 
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Fig. 14. Coherence (a) and Shannon entropy (b) of the asymptotic state in the vicinity 
of a singlet-doublet crossing for different temperatures as given in the legend. The other 
parameter values are D — A, fi = 0.982 tJo, and 7 = lO^^cJo. 




Fig. 15. Coherence of the asymp- 
totic state in the vicinity of a sin- 
glet-doublet crossing for F — 0.013 
(a) and F = 0.015029 (b): exact 
calculation (full line) compared to 
the result of a three-level descrip- 
tion (dashed) of the dissipative dy- 
namics. The other parameter val- 
ues are D = 4, f2 = 0.982 wo, and 
7 = IQ-'^uJo- 
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temperatures fceT <C 2b and much less than unity for fceT 3> 2b (Figs. 14a, 15). 
This means that the crossing of the chaotic singlet with the regular doublet leads 
to an improvement of coherence if the temperature is below the splitting of the 
avoided crossing, and to a loss for temperatures above the splitting. This phe- 
nomenon amounts to a chaos-induced coherence or incoherence, respectively. The 
corresponding Shannon entropy (Fig. 14b) assumes approximately the value Inn 
for n incoherently populated states. Thus outside the crossing, we have S w ln2 
and at the center of the crossing the entropy exhibits a significant temperature 
dependence. 

The crucial role of the decay via states not involved in the three- level crossing 
can be demonstrated by comparing it to the dissipative dynamics including only 
these three levels (plus the bath). At the crossing, the three-state model results 
in a completely different type of asymptotic state (Fig. 15). The failure of the 
three-state model in the presence of dissipation clearly indicates that in the 
vicinity of the singlet-doublet crossing, it is important to take a large set of 
levels into account. 
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5 Signatures of chaos in the asymptotic state 

Phase-space representations of quantum mechanics, hke the Husimi or the Wig- 
ner distributions, help to reveal the structures of the corresponding classical 
phase space [7,58-62]. In particular, for the case of regular classical dynamics, 
the Husimi function of a (quasi)energy eigenstate is localized along the corre- 
sponding quantizing torus; for chaotic motion, it is spread over the entire chaotic 
layer. If the classical dynamics is mixed, quantum-mechanical states can be clas- 
sified as regular or chaotic according to their distribution in phase space [62]. 
Moreover, the phase-space representation of the asymptotic state of a dissipative 
quantum map exhibits the structures of the corresponding classical attractor [23] . 
However, these analogies find their limit in the Heisenberg uncertainty princi- 
ple. It does not allow for arbitrarily fine phase-space structures and results in 
smearing on action scales below h. 

5.1 Classical attractor 

To describe the classical dissipative dynamics of the driven double well, we add 
an Ohmic friction force F^ = —jp to the conservative equations (4), (5), 

X = -p, (59) 

m 

P^-IP Q^- (60) 

Friction destroys the time-reversal symmetry (21) of the conservative system. 
Accordingly, dissipation breaks the reflection symmetry of the phase-space por- 
trait with respect to the x-axis, found at zero phase of the driving (cf. Fig. 6). 

A constituent feature of dissipative flows is the net exponential contraction 
of phase-space volume. Therefore, the dynamics is asymptotically confined to 
an attractor, a measure-zero manifold in phase space to which all trajectories 
starting from within the surrounding basin of attraction converge. For periodi- 
cally driven dissipative systems, the attractor is in general also time-dependent 
with the period of the driving and is adequately described by a stroboscopic 
map [63-65]. 

Depending on the parameters that control the dissipative flow, an attractor 
can consist of fixed points, limit cycles, or manifolds of fractal dimension, less 
than that of phase space. An adequate concept to characterize the geometry of 
an attractor is the Hausdorff dimension dn defined, for example, in Ref. [3]. It 
typically increases with decreasing contraction rate, so that strange attractors 
are expected to occur in the regime of weak dissipation of a system that in 
absence of friction, is chaotic. 

The Hausdorff dimension of the classical attractor for the driven double well 
with dissipation, Eqs. (59), (60), at the parameter values F — 0.09 and fi = 
0.9 Wo, is shown in Fig. 17 for various values of the friction 7. Although the 
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Fig. 16. Stroboscopic classical phase 
space portrait at f = nP, of the 
dissipative harmonicaUy driven quar- 
tic double well, Eqs. (59), (60), for 
the driving amplitude F — 0.09 and 
frequency J7 = 0.9 uiq. The friction 
strength is 7 = 0.3ci;o (a), 0.2 cjo (b), 
0.03 ojo (c). In panels (a) and (b) the 
stroboscopic portrait is marked by a 
full dot and the broken lines show the 
corresponding limit cycles. 




Fig. 17. Hausdorff dimension of the 
classical attractor of the dissipative 
harmonically driven quartic double 
well, Eqs. (59), (60), for F = 0.09, 
n = 0.9uJo- 



attractor itself is periodically time dependent with the period of the driving, its 
Hausdorff dimension dfj remains nearly constant [63]. Near 7 « 0.06 wq, with 
decreasing 7, the classical dynamics undergoes a transition from regular motion 
to chaos, manifest in a corresponding transition from limit cycles (Fig. 16a, b) to 
a strange attractor (Fig. 16c) and a concomitant jump in d^i- At the same values 
of F and f2, the regular islands near the potential minima of the corresponding 
undamped dynamics have already completely dissolved in the chaotic sea. 
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Fig. 18. Contour plot of the Husimi function of the quantum attractor (full lines) 
at i = nP, n -^ cx), superposed on the corresponding classical phase-space portrait, 
Fig. 16. The parameter values F = 0.09, fi — O.Qojq, 7 ~ 0.3ujo (a,d), 0.2ciJo (b,e), 
0.03 uo (c,f) are as in Fig. 16. The effective action is D = 6 (a-c) and D = 12 (d-f). 
The rectangle in the lower left corner has the size of the effective quantum of action 
Reff = 1/8D. 
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Fig. 19. Wehrl entropy of the 
asymptotic state of the dissipative 
quantum map for different values 
of the effective quantum of action 
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5.2 Quantum attractor 

The "quantum attractors" , i.e., the asymptotic states of the dissipative quantum 
dynamics, for example in a Husimi representation, resemble the correspond- 
ing classical attractors up to coarse graining (Fig. 18). Correspondingly, the 
qualitative transformation from limit cycles to a strange attractor is visible 
in the asymptotic quantum distribution, but proceeds continuously. Although 
the asymptotic state for 7 — 0.2 wq (Fig. 18b, e), is still concentrated near the 
fixed points of the classical stroboscopic map, it covers a broader phase-space 
area that already anticipates the shape of the strange attractor. Reducing the 
effective quantum of action h^g — 1/8D allows for a sharper resolution of the 
underlying classical structures in the Husimi functions, as expected. 

Like the attractors of the dissipative classical dynamics (Fig. 16), their quan- 
tum-mechanical counterparts have lost the reflection symmetry with respect to 
the X-axis, in contrast to the Husimi representations of the Floquet states in 
absence of dissipation (cf. Fig. 7). This symmetry breaking is caused by finite 
off-diagonal elements of the asymptotic density matrix in Floquet representa- 
tion, since diagonal representations share the symmetries of the basis. Thus, 
off-diagonal matrix elements play a significant role for the asymptotic state. 
This demonstrates that a description within a full rotating-wave approximation 
is insufficient, since it would result in an asymptotic state diagonal in the Floquet 
representation [51-53]. 

Because the quantum attractor, in contrast to the classical one, has support 
all over phase space (or a region of finite measure), we cannot characterize it 
by a Hausdorff dimension. A more suitable measure for the extension of the 
quantum attractor is a phase-space version of the Shannon entropy, the Wehrl 
entropy [62,66,67] 



SQ 



of its Husimi representation 



dxdp 
2Trh 



Q{x,p)\n[Q{x,p)] 



(61) 



Q(x,p) = {x,p\ps\x,p), 



(62) 
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where \x,p) denotes a coherent state centered at (x,p) in phase space. Its expo- 
nential, exp(S''5), gives approximately the number of coherent states covered by 
the Husimi function. Thus, the occupied phase-space area is 27r?iexp(S''5). The 
Wehrl entropy of the asymptotic state for our numerical example is depicted in 
Fig. 19 for various values of ?ioff = 1/813. It grows with decreasing friction 7, 
reflecting the increasing dispersion of the Husimi functions. In the semiclassical 
regime, i.e., for a sufficiently large value of the effective action 13, we observe a 
kink-like behavior of the entropy near 7 « 0.06 wq, where the classical attractor 
undergoes the transition mentioned above, from a set of isolated fixed points to 
a strange attractor. 

Note that for 7 <; 0.1 wq, the Markov approximation becomes inaccurate, 
since 7 is then of the order of the mean level spacing and the condition (43) is 
violated for at least part of the transitions between Floquet states. Nevertheless, 
we obtain the qualitative behavior which we expected from classical considera- 
tions. 



6 Conclusion 

For the generic situation of the dissipative quantum dynamics of a particle in a 
driven double-well potential, classical chaos plays a significant role for the co- 
herent dynamics. Even for arbitrarily small driving amplitude, the separatrix is 
replaced by a chaotic layer, while the motion near the bottom of the wells re- 
mains regular. Nevertheless, the influence of states located in the chaotic region 
alters the splittings of the regular doublets and thus the tunnel rates, which is 
the essence of chaotic tunneling. We have studied chaotic tunneling in the vicin- 
ity of crossings of chaotic singlets with tunnel doublets under the influence of 
an environment. As a simple intuitive model to compare against, we have con- 
structed a three-state system which in the case of vanishing dissipation, provides 
a faithful description of an isolated singlet-doublet crossing. Dissipation intro- 
duces new time scales to the system: one for the loss of coherence and a second 
one for the relaxation to an asymptotic state. Well outside the crossing, both 
time-scales are of the same order, reflecting an effective two-state behavior. The 
center of the crossing is characterized by a strong mixing of the chaotic state 
with one state of the tunnel doublet. The high mean energy of the chaotic state 
introduces additional decay channels to states outside the three-state system. 
Thus, decoherence becomes far more effective and, accordingly, tunneling fades 
out much faster. 

The study of the asymptotic state, the quantum attractor, demonstrates 
clearly that a three-state model of the singlct-doublct crossing is insufficient 
once dissipation is effective. This is so because the coupling to the heat bath 
enables processes of decay and thermal activation that connect the states in the 
crossing with other, "external" states of the central system. In the presence of 
driving, the asymptotic state is no longer literally a state of equilibrium. Rather, 
incoherent processes create a steady flow of probability involving states within 
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as well as outside the crossing. As a result, the composition of the asymptotic 
state, expressed for example by its coherence trp^, is markedly different at the 
center of the crossing as compared to the asymptotic state far away from the 
crossing, even if that is barely visible in the corresponding classical phase-space 
structure. 

With increasing driving amplitude, in absence of dissipation, even the dynam- 
ics near the bottom of the wells becomes fully chaotic. This has striking con- 
sequences for the corresponding dissipative classical dynamics: For sufficiently 
weak dissipation, it remains chaotic, but for strong friction it becomes regular. 
Accordingly, the geometry of the classical attractor is fractal or regular, re- 
spectively. We have observed the signatures of this qualitative difference in the 
asymptotic state of the corresponding quantum dynamics. However, in contrast 
to the sudden change of the classical behavior, the quantum attractor under- 
goes a smooth transition: The structure of the strange attractor is already felt 
by the Husimi function for parameter values where the classical attractor con- 
sists only of two isolated fixed points. For the observation of these semiclassical 
structures, off-diagonal matrix elements of the asymptotic state in Floquet basis 
proved crucial. This clearly indicates that a full rotating-wave approximation 
must fail. 

Many more phenomena at the overlap of chaos, tunneling, and dissipation 
await being unraveled. They include four-state crossings formed when two dou- 
blets intersect, chaotic Bloch tunneling along extended potentials with a large 
number of unit cells instead of just two, and the influence of decoherence on a 
multi-step mechanism of chaotic tunneling. These phenomena are typically ob- 
served in the far semiclassical regime, which requires to take very many levels 
into account. A semiclassical description of the dissipative quantum system may 
circumvent this problem. 
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